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Abstract 

We describe the range of the spherical Radon transform which 
evaluates integrals of a function in IR" over all spheres centered on a 
given sphere. Such transform attracts much attention due to its ap- 
plications in approximation theory and (thermo- and photoacoustic) 
tomography. Range descriptions for this transform have been obtained 
recently. They include two types of conditions: orthogonality condi- 
tion and, for even n, a moment condition. However, it was found out 
later that, in any dimensions, the moment condition follows from the 
orthogonality one, and therefore can be dropped. In terms of Darboux 
equation, which describes spherical means, it indirectly implies that 
solutions of certain boundary value problems in a domain automati- 
cally extend outside of the domain. In this article, we present a direct 
proof of this global extendibility phenomenon for Darboux equation. 
Correspondingly, it delivers an alternative proof of the range charac- 
terization theorem. 



1 Introduction 

Let / be a function defined on IR". Consider the spherical mean Radon 
transform 

nf{x,t):=— I f{x + te)dA{0), 

J 

where dA[9) is the area measure on the unit sphere 5""^ C IR" and Un is 
the total measure of the unit sphere. This transform has been intensively 
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investigated in the literature due to its applications in PDEs, medical imaging 
and geophysics (e.g., [2, 6-9, 14, 15]). 

It is well known that G{x,t) = TZf{x,t) satisfies the Darboux equation 
[8, 15] 

r VG{x,t):=0, e IR'' X 1R+, , . 

\ G(x,0) = /(x), Gi(x,0) =0. 

Here 1R+ := (0, oo) and V is the Darboux operator 

77 — 1 

V := + —^dt - A.. (2) 

Moreover, Aisgeirsson's theorem [8, Ch.6] states that any global solution 
G{x,t) of (1) is the spherical means of its initial value: G{x,t) = {TZf){x,t). 

The problem of recovering a function f of n variables from the function 
TZf{x,t), which depends on n + 1 parameters, is obviously overdetermined. 
Thus, it is natural to consider the restriction of TZf{x, t) to an n-dimensional 
surface in the (x, t)-space IR*^ x 1R_|_. Different problems in analysis are related 
to different types of such surfaces. 

Problems concerning the spherical mean operator restricted to cylindrical 
surfaces in IR" x 1R_|_ have attracted a lot of attention recently. By cylindrical 
surfaces we mean those of the form T := S x 1R_|_ where S* is a hypersurface in 
IR" . The problem of recovering finctions from their spherical means restricted 
to such surfaces has interesting applications in analysis, in particular, in 
approximation theory (see [4,17]). But perhaps more importantly, it serves 
as a mathematical model for thermoacoustic tomography, a novel medical 
imaging method (we refer the readers to [2, 12, 16] for detailed explanation). 
In that application, the restriction g = TZsif) of T^if) to the set T is the 
measured data, while / is the image to be determined. 

In this article, we are concerned with the problem of characterizing func- 
tions g that belong to the range of the operator TZs- 

Let us describe the problem in more precise terms (for the detailed ex- 
position, we refer the reader to the articles [3,5,13]). First of all, in our 
considerations, the hypersurface S will be the unit sphere centered at the 
origin. We assume that / e C°°(IR") is supported inside the closure B of 
the unit ball i? = {|x| < 1} (the class of such functions / will be denoted by 

G^m- 

We consider the problem of characterizing all functions g{x,t) on the 
cylinder F such that g = TZsf for some function / G G^{B). 

Some necessary conditions are almost obvious. First, the function g{x,t) = 
TZsf{x,t) must be smooth on F. Second, g{x,t) must vanish for all t > 2 
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and also vanish to infinite order at t = 0. Thus, g must satisfy smoothness 
and support conditions: g G C^{S x [0,2]). 

Another, less trivial, necessary condition can be derived from characteri- 
zation of the spherical means by Darboux equation (1). This condition, which 
we call orthogonality condition (see Theorem 2.3), is of Fredholm alternative 
type. It follows from the existence of a solution of Darboux equation inside 
the cylinder F (the fact that is a sphere is not important here). 

The first complete range description for n = 2 was obtained in [5]. It was 
proved that a function g belongs to the image of C^{B) under the operator 
TZs if and only if it satisfies, besides the smoothness and support conditions 
and the above orthogonality condition, an additional moment condition. Ear- 
lier, the necessity of the moment condition was observed in [18]. 

Further step in higher dimension was taken in [13]. There a range charac- 
terization was obtained for odd dimension and for a transform related to the 
wave equation, rather than Darboux equation. The conditions in [13] did not 
involve the moment conditions. A complete description of the range for the 
spherical mean transform in any dimension was obtained in [3]. As in [5], the 
necessary and sufficient conditions in [3] fall into two groups: orthogonality 
and moment conditions, although the moment conditions were needed for 
even dimensions only. 

It had remained unclear whether the moment condition is really needed 
for even n, till the recent article [1]. It was proved there that, regardless 
of parity of the dimension, the moment conditions follow from smoothness, 
support, and orthogonality conditions and therefore can be dropped. 

The above result can be immediately translated to the language of Dar- 
boux equation. Namely, on one hand, it was proved in [3] that the orthog- 
onality condition for the data g{x., t) is in fact the condition of existence of 
a solution G = of Darboux equation (1) inside the cylinder F with the 
boundary data G^{x, t) = g{x, t) on F and proper decay when t —>■ +oo. On 
the other hand, spherical means are global solutions of Darboux equation. 
Therefore, g is in the range of transform TZs means that G"*" extends as a 
global solution to the entire space IR"^. 

The possibility of such an extension seems to be an interesting phe- 
nomenon by itself and one may wish to have its direct proof. In fact, such 
a proof was found in [3], but in odd dimensions. In this article, we mod- 
ify the construction of [3] to extend it to all dimensions. Correspondingly, 
we obtain an alternative proof, universal for all dimensions, of the range 
characterization theorem from [1]. 
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2 Main result 



Let us recall here that B is the unit ball centered at the origin, S = dB is the 
unit sphere, and V = S x IR+. We also denote by the Bessel function of 
the firts kind of order ^ and jfi{u) = u~^J^{u) the corresponding normalized 
Bessel function. The notation C^{VL) will stand for smooth functions in IR" 
with the support in f2. 

2.1 Formulation of main result 

The goal of this article is to present a direct proof of the following result 
from [1]: 

Theorem 2.1. Let g he a function defined on the cylinder T . Then there 
exists f G C^{B) such that g = TZs{f) if and only if the following conditions 
hold: 

a) Smoothness and support conditions: g G C^{T) and g{x,t) = 
when t > 2. 

b) Orthogonality condition: Let —}? he an eigenvalue of the Dirichlet 
Laplacian on B and tpx o corresponding eigenfunction. Then 

J g{x,t)d^^<^x{x)jn^{Xt)t''-^dtdA{x) = 0, 
r 

where Vx is the outward normal to S at x. 

Remark: Since S is the unit sphere, the Dirichlet eigenfunctions ip\ can 
be written in polar coordinates as follows: 

•^xirO) = j n+n^-2 (\r)Yrr,(9). 

where Ym is a spherical harmonics of degree m. The Dirichlet condition for 
(fx on the unit sphere S requires that j(„+.m_2)/2(A) = 0. Choosing Ym = 
Ym, k, k = 1, ■ ■ ■ , d{m), elements of the basis in the space of all harmonics 
of degree m, one can write the condition b) in the equivalent form: 

b') J g{e,t)jn^{Xt)t''-^Y^,k{0)dtdA{e) =0. 
r ^ 

This can be rephrased as follows: 
b") ?m,fc(A) = 0, for all zeros A of the Bessel function j n+m-2 (A), where 
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is Fourier-Bessel transform of gjn,k 
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Theorem 2.1 can be reformulated in terms of Darboux equation. Namely, 
since the spherical means G = TZf is the unique solution for the Darboux 
equation (1), Theorem 2.1 is equivalent to: 

Theorem 2.2. Let g be a function defined on the cylinder T. Then the 
following statements are equivalent: 

i) There exists f G C^{B) such that the following problem has a solution: 



a) The conditions a) and b) of Theorem 2.1 hold. 

The proof of the implication i) =^ ii) is quite simple and can be found in 
[3]. In the rest of this article, we will prove the converse implication. 

2.2 Orthogonality condition and existence of the in- 
ternal solution 

In this subsection, we recall a result from [3], which is our starting point for 
the proof of Theorem 2.2. This result interprets the orthogonality condition 
as the existence of some internal solution for the Darboux equation as follows: 

Theorem 2.3. Let T > 0. Consider the following backward boundary initial 
value problem for Darboux equation: 



where Q C IR" is a bounded domain with smooth boundary and g G G^{dQ x 
[0,r]). Then the following statements are equivalent: 




(3) 
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a) The solution G G C°°{Q x IR+) is smooth at t = and G{x,t) = for 
t > T. 

b) The boundary data g satisfies the orthogonality condition 

g{x,t)d^^(pxix)jn_2{Xt)e-^dA{x)dt = 0, (4) 
an 

for all pairs of Dirichlet eigenvalue-eigenfunctions (A, ipx) ■ 

The detailed proof can be found in [3], and here we only briefly explain 
its main idea. The implication a) =^ b) immediately follows from Stokes 
formula applied in x variable and integration by parts with respect to t. Vice 
versa, assume that the orthogonality condition (4) holds. The equation is 
nonsingular for all t > 0, and the unique solution exists in i? x (0,T]. One 
needs only to show that G is not singular at t = and Gt{x, 0) = 0. This is 
done by applying Stokes formula in the domain B x [e,T], e > 0, using the 
orthogonality condition (4) and letting e 0. . 

Remark: The above orthogonality condition can be understood as a 
Fredholm alternative statement for solvability of the boundary value problem 
for the equation VG = in the solid cylinder Q x 1R+. Indeed, denote 
by Pg{x,t) the harmonic extension (Poisson integral) in x of the boundary 
data g{x,t),x e S. By setting H = G — Pg, one can rewrite the boundary 
value problem 1 for G as a boundary value problem for the nonhomogeneous 
equation 

71 — 1 

VH = -{d^ + —^dt)Pg, 

with zero boundary data H{x,t) = 0, x ^ S. Now, Fredholm alterna- 
tive claims that for solvability the right hand side of the equation must 
be orthogonal to all solutions of the homogeneous adjoint equation. Since 
Darboux operator with Dirichlet boundary conditions is self-adjoint, it suf- 
fices to check the orthogonality of the right hand side to separable solutions 
ux{x,t) = jx{t)(px{x), which constitute a complete system of solutions. Then 
Fredholm alternative condition becomes 

uxix,t)Pg{x,t)t''~^dtdx = 0, 

which is just our orthogonality condition (4), if uses the Stokes formula to 
replace the integration in x over Q by the surface integration over dQ. 
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2.3 Theorem 2.2 as extendibility theorem for Darboux 
equation 

Let C = B X IR+,r = S X IR_|_, and suppose that the conditions a) and b) 
of Theorem 2.1 are staisfied. Consider the backward initial boundary value 
Darboux problem 

VG{x,t) = 0, {x,t) G C, 
G{x, t) = g{x, t), {x, t) G r, 
Glx,t) = 0, {x,t) eBx [2,oo). 

Theorem 2.3 for the case Q = B and T = 2 shows that the solution G 
C°°(C) of this problem is non-singular at t = and Gt{x, 0) = 0. On the other 
hand, if a global solution G for (3) exists, it also solves Darboux equation 
in the cylinder C and satisfies the same boundary and initial conditions: 
G{x, t) = g{x, t), {x, t) G r, and G{x, t) = for t > 2. Due to the uniqueness 
of the solution, G and G^ coincide on C. That is, the existence of G is 
equivalent to the global extendibility for G^. 

The above argument shows that the implication ii) =^ i) in Theorem 2.2 
can be reduced to the following extendibility result: 

Theorem 2.4. Let G"*" G C°°{C) be a solution of Darboux equation 

VG+{x,t) = 0, {x,t) G C, 

such that Gt{x,0) = and G~^{x,t) = for t > 2 and x E B. If the 

boundary value g = G'^\r belongs to G^{S x [0,2]), then G^ extends to 

a global solution G of the Darboux equation (3), with the C°° initial value 
f*{x) that is G~^{x,0) extended by zero outside B. 



3 Proof of Theorem 2.4 

The proof rests on the following two propositions. The first one says that if 
the initial data / of the internal solution G^ vanishes to infinite order on the 
unit sphere then the extendibility holds. 

Proposition 3.1. [3] Let G'^ G C°°(C) be a solution of Darboux equation 

VG^{x,t) = 0, {x,t)eC 

and G^{x,t) = for all t > 2 and x E B. If f{x) := G{x,Q) vanishes to 
infinite order on the unit sphere \x\ = 1, then G^ extends to a global solution 
G{x,t) of Darboux equation. This global solution is given by the spherical 
means G{x,t) = {TZf*){x,t) of the function f* in IR" obtained from f by the 
zero extension outside of the unit ball. 
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The second proposition proves that the smoothness and support condi- 
tions for the boundary data g imply required vanishing to infinite order of 
the function / on the sphere S. 

Proposition 3.2. Let G C°°(C) be a solution of Darboux equation 

VG{x,t) = 0, V(2;,t) e C 

such that G{x,t) = for t > 2 and x & B. If the boundary value g = G^\r 
belongs to G^{S x [0,2]) then f{x) = G{x,0) vanishes to infinite order on 
the sphere \x\ = 1. 

We will prove Propositions 3.1 and 3.2 in separate sections. Meanwhile, 
we derive Theorem 2.1 from these two propositions. 

The proof of Theorem 2.4 is just simple combination of Propositions 
3.1 and 3.2. Indeed, let G^{x,t) be the internal solution as in Theorem 2.4. 
Then Proposition 3.2 implies that f{x) := G~^{x, 0) has zero of infinite order 
on the unit sphere = 1. By Proposition 3.1, the (globally defined) spher- 
ical means G{x,t) := TZf*{x,t) coincide with G~^{x,t) inside the cylinder F 
and therefore provide the extension of the internal solution G~^ . This proves 
Theorem 2.4 and hence proves the equivalent Theorem 2.1. 

4 Proofs of the propositions 
4.1 Proposition 3.1 

As we have already mentioned. Proposition 3.1 is proven in [3]. We will 
present the proof here for the sake of completeness. 

Since f{x) = G^{x, 0) vanishes to infinite order on the boundary of B, the 
zero extended function /* belongs to C°°(1R"). Then the natural candidate 
for the extended solution is given by the spherical means of /*: 

G = nf*. 

This function is globally defined and, since /* is smooth, belongs to C°°(lR"x 
IR+). It is a global solution of Darboux equation and our goal is to prove 
that G'^ and G coincide in the solid cylinder C := B x [0, oo). 

First observe that both solutions share the same initial data at t = 2: 

G(x,2) = G+(a;,2) = Gt{x,2) = Gt{x,2) = 0. 

Then, due to domain of dependence theorem [8, p. 696], 

G{x, t) = G+{x, t) = for all {x, t) G /C+, (5) 
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where K,^ is the upward characteristic cone 

;C+ := {(x, t) G X [0, 2] : f - |x| > 1}. 

Moreover, G{x,t) and G~^{x,t) also share the initial data at t = 0: 

G{x,0) = G+{x,0) = f{x),Gt{x,0) = Gi{x,0) = 0. 

Therefore, again by the dependence domain theorem, they coincide in the 
downward characteristic cone with the base B x {0}, 

/C" := {{x,t) eBx [0, 1] : |x| + t < 1}. 

Hence, the difference U{x,t) := G{x,t) — G~^{x,t) vanishes in the union 
/C = /C"*" U /C~ of the two cones. Besides, since both G and G~^ vanish for 
t > 2, so does U. 

Since U satisfies Darboux equation inside the cylinder C, its Fourier- Bessel 
transform 

oo 

U{x,a) = j U{x,t)jn^{at)e-^dt 



satisfies Helmholtz equation: 

AxU{x,a) = —a'^U{x,a), Vx G B. 
Hence, U{x, a) is real analytic with respect to x E B. Also 

U{x,t) = 0,x G /C 

and U{x,t) = for t > 2. The union /C U (S x [2,oo)) contains the entire 
ray {(0, t) : < t < oo} and hence after taking Fourier-Bessel transform one 
conludes that U{0,a) = 0. Since U is smooth, the same argument can be 
applied to D^U to obtain: 

D'^U{0,a) = D^{0,a) = 0. 

Thus, U{., X) vanishes to infinite order at x = 0. Since U{., X) is real- analytic, 
one concludes that U{x,a) = 0, x E B for all a. Taking inverse Fourier-Bessel 
transform, we obtain [7 = and therefore G = G^ in C. This completes the 
proof of Proposition 3.1. 
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4.2 Proof of Proposition 3.2 

We want to prove that / := G^(.,0) vanishes on the sphere S to infinite 
order: 

= 0, \x\ = 1. 

First of all, recall that G{x, 2) = Gt{x, 2) = implies = in the upward 
characteristic cone }C~^: 

G^{x,t) =0, |x| < t - 1, 1 < t < 2. 

In particular, G~^ vanishes at the vertex of the cone /C"*": 

G+(0,1) = 0. 

Since G"*" is smooth in the solid cylinder C := B x [0, oo), the same conclusion 
holds for for all derivatives of G"*": 

D|DfG+(0,l) = 0. (6) 

Here j = 0, 1, ... and (] is an arbitrary multiindex. 

Now we relate G"*" to the spherical means of the initial value f{x). Albeit 
we cannot assert so far that G~^ = TZf (which is in fact our final goal), we can 
claim that the two functions coincide at least in the downward characteristic 
cone /C~: 

/C" = {\x\ <l-t, < t < 1}. 

Indeed, both G^ and TZf solve the equation VG = on C and share the same 
initial values on B: G{x,0) = f,Gt{x,0) = 0. The conclusion now follows 
from the domain of dependence argument. 

Since (0, 1) is the vertex of /C~, due to (6), we conclude that 

(D^Df7^/)(0,l) = 0. (7) 

Now our aim is to derive from (7) that f{x) vanishes on the sphere \x\ = 1 
along with all derivatives. 

To this end, first observe that condition (7) is invariant with respect 
to action of the orthogonal group 0{n) and hence it holds for any term 
fm,k{^)r''^ym,k{0) in Fourier series of /. 

oo d(m) 

fix) = fire) = J2Y. fmAry'^ymAO). (8) 

m=0 k=l 

where r = |x|, |^| = 1 and Y^^kyk = l,...,(i(m) is the orthonormal basic in 
the space of all spherical harmonics of degree m. Due to smoothness of / 
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in the closed ball B , the series (8) converges uniformly with all derivatives 
and hence it suffices to prove that each term vanishes on the unit sphere to 
infinite order. 

Thus, we can assume that / is just a single term: 

fix) = fm{r)Pm{x), (9) 

where 

PUx) = r'"F^,fc(e) 

is a spatial harmonic of degree m. To prove that / vanishes to infinite order 
for \x\ = 1, it suffices to prove that all the derivatives fm vanihs at t = 1: 

/(^■)(i) = o,j = o,i,... 

We will prove this by constructing a system of linear equations that these 
numbers satisfy. 

Lemma 4.1. The following identities hold 
1. For any i > 0, 

{^.Cmfmm = ^, (10) 

where Cm is the following differential operator of order m: 

^-=\{i ^^,}_. 4+i)- (11) 



s=l 



2. For any I > 0, 

(Q,n)'(/m)(l)=0, (12) 

where Qm is the following differential operator of order 2: 

Qm = dl + dr. (13) 



Proof: It will be convenient to introduce the operator 

dim) 



{9){r) = Y.gm,k{r)r^Y^{e). 



k=l 
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that projects onto spherical harmonics of degree m. We observe that differ- 
entiation in X and the transform TZ commute: 

Df7^/(x,^) = 7^(D^/)(x,^). 

Since the spherical mean with the center at a; = and radius t = 1 is exactly 
the projection onto order zero harmonics (constants), identity (7) now reads 
as 

-^7ro(D^/)(l) = 

for all i G and all multiindices f3. 

The projection tcq of the derivatives of / of the form (9) was computed 
in [9, formula (2.10)]: 

MD^f) = [D^Prr) {Crrdm) , m = = A + ■ ■ ■ + 

where the differential operator Cm is defined in (11). Since Pm is a polynomial 
of degree m, we can choose in this formula the multiindex (3 so that D^Pm 
is a non-zero constant. By allowing the index i to be arbitrary, we arrive at 
the identity (10). 

As for the identity (12), it comes from the equation T>G^ = 0, which 
means 

BG+{x,t) = A^G+{x,t), \/{x,t) G C, 
where B is the Bessel operator acting on t— variable 

Iterating the above identity, one obtains 

B^G+{x, t) = A^G+{x, t), V(x, t) G C. (14) 

Since G~^ G C°°{C), the above equality holds up to the boundary T. In 
particular, since G{x,t) = g{x,t) for all {x,t) G F and G{x,0) = f{x), we 
have 

B^g{x,0) = A7(x), Vx G S. (15) 
Since g vanishes to infinite order at t = 0, one concludes that for all x G S 

A7(a;) = 0,/ = 0,1,... 

Now, taking into account that the Laplacian A acts on m*^-harmonic term 
as the operator 

o n — 1 + 2m „ 

Qm := + dr, 

r 
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we arrive at 

QUm = 0,1 = 0,1,... (16) 

■ 

Let us apply Lemma 4.1 when indices j and / run independently from 
to m — 1. We can write the differential operators in (10) and (12) in the form 

-f-£rn = A,oir) + AiUr)-^ H h Aii+mir) ^ 



ar^ ar ar^^^ 

and 

ar ar '' 

Consider the vector 

i^:=(/™(l),/;(l),---,/i^-^Hl)), 

which consist of the first 2m derivatives (including that of order 0) of fm at 
the point r = 1. 

Let i = 0,1, ■ ■ ■ ,m — l and / = 0, 1, ■ ■ ■ , m — 1 in Lemma 4.1, we conlcude 
that F satisfies the following 2m x 2m linear system: 

^0,0-^0 + ^0,1-^1 + 1- ^0,m.-^m = 

^1,1-^1 + ^1,2-^2 + + ^l,m+l-^m+l = 



^m— l,m — 1 ~l~ — l,m-^m ~l~ ~l~ ^m— 1,2m— l-f2m— 1 

{ Bo,oFo = (17) 
Bi^iFi + ^1,2^^2 = 

B20F2 + ^2,3-^3 + B2AF4 = 



Bm-l,m~lFm-l + + i?m-l,2m-2-^2m-2 = 0. 

Here the matrix coefficients are Aij = Aij{l), Bij = Bi j{l). 

Lemma 4.2. The linear 2m x 2m— system (17) is nondegenerate. 

We will prove this lemma later. Assuming that the lemma is proven, we 
can complete the proof of Proposition 3.2. Since the system (17) is nonde- 
generate, one concludes that the first 2m derivatives of fm vanish: 

fli\l) = F^ = 0,0<3<2m-l. 

To obtain the vanishing of higher order derivatives, we will exploit higher 
values for the index i in (10). Choosing i = m in (10) results in shifting of 
the vector of the unknowns to the right: 

+ Am,2in-lF2m-l + Am,mF2m = 
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which along with Fq = ■ ■ ■ = F2m-i = imphes = because Am,m 7^ 0. 
Then the next choice i = m + 1 leads to -F2m+i = 0. Proceeding this way by 
taking successively j = m, m + 1, m + 2, ■ ■ ■ one obtains = for all u > 0. 
This completes the proof of Proposition 3.2. ■ 

Remark: As it was mentioned earlier, Proposition 3.2 was proved for 
odd n in [3]. The proof used Weyl and Poisson-Sonine integral transforms, 
applied to the solution G^{x, t) in t and x variables correspondingly. To have 
control over the derivatives of G^{x, 0) on the sphere |x| = 1 one needs the 
inverse transforms to be local (differential) operators which is the case only 
in odd dimensions. That is why the proof in [3] did not generalize to even 
dimension. 

Proof of Lemma 4.2 

Let us denote by Ai the ith row of the matrix of the first m equations 
in the system (17) and by Bi the ith row from the second group of m linear 
equations in (17). 

One observes that the vectors Ai,i = 0, ■ ■ ■ ,m — 1 are linearly indepen- 
dent, as Aij = for all j > m + i and Ai^rn+i 7^ 0. 

Now we will use induction with respect to the length of the system of 
the vectors. Namely, we will show that on each step successive addition of 
vectors Bq,--- , Bm~i to the set {A^ , A^"^} does not violate the linear 
independence of the set obtained on the previous step. Then, in m steps, we 
will obtain the linear independence of the entire system. 

Thus, our inductive assumption is that the system 

Sp := {Ai, ■ ■ ■ , Am-i, Bq, ■ ■ ■ , Bp^i} 

is linearly independent for some p < m — 1. Now we want to check that 
it remains linearly independent after adding the next vector Bp. In other 
words, the vector Bp is linearly independent from the set Sp. 

To this end, it suffices to find a vector Vp E IR^™ that is orthogonal to Sp 
but not to Bp-. 

{Ai,Vp) = 0,z = 0, ■ ■ ■ ,m - 1, 
{Bj,Vp) = 0,j = O,--- ,p- 1, 
{Bp,vp) + 0. 

Indeed, we take the function '^p[f) = and construct the vector of 

successive derivatives at r = 1: 

vp = {%{!),■ ■■ 
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Recall that 



n 



1 



n + 2(m — s) dr 



-r— + 1 . 



Since all the first order differential operators in the above product commute, 
we can rewrite 



n 



n + 2(m — s) dr 
A simple observation gives 



d ^ 
r— + 1 



1 d 
-r 



n + 2p dr 



^ -r^ + 1 1 ^Jr) = 0,Vr > 0. 



n + 2p dr 



Therefore, 



£^*p(r) = 0, Vr > 0. 



This, in particular, implies 

(J-Cm'^p^ (1) = 0, ^ = 0,1,--- ,m-l. 

By the definition of Ai, this is exactly equivalent to the first group of equa- 
tions in (17): 

{Ai,Vp) = 0,z = 0, ■ ■ ■ ,m - 1. 

Here, as we defined, Vp = (\E'p(l),--- is the vector of successive 

derivatives of evaluated at r = 1. 

The second group of the equations in (17) comes, by the construction, 
from the iteration of the differential operator 

r 

evaluated at r = 1. Straightforward computation yields 

' C';r-"-2(P+'), < / < m - 1 



2 , n-l+ 2m 
ot H Or 



^>p{r) 



0, m — p<l<m — 1. 



Here Ci are nonzero constants. 

Substituting r = 1 and recollecting the original construction of the vectors 
Bq, - ■ ■ , Bm-i one is led to 



{Bo,Vp) 



{Bp^i, Vp) = 0, {Bp, Vp) = Cm-i-p 7^ 0. 



This completes the proof of Lemma 4.2 and thus finishes the proof of Theorem 
2.1. ■ 
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